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ABSTRACT 

NGC 3576 (G291.28-0.71; I = 291.3°, b = -0.7°) and NGC 3603 (G291.58-0.43; I = 
291.6°, b = -0.5°) are optically visible, luminous HII regions located at distances of 3.0 
kpc and 6.1 kpc, respectively. We present 3.4 cm Australian Telescope Compact Array 
(ATCA) observations of these two sources in the continuum and the H90a, He90a, 
C90a and HI 13/3 recombination lines with an angular resolution of 7" and a velocity 
resolution of 2.6 km s -1 . All four recombination lines are detected in the integrated 
profiles of the two sources. Broad radio recombination lines are detected in both NGC 
3576 (AV FWHM > 50 km s" 1 ) and NGC 3603 (AVfwhm > 70 km s" 1 ). In NGC 
3576 a prominent N-S velocity gradient (~30 km s _1 pc -1 ) is observed, and a clear 
temperature gradient (6000 K to 8000 K) is found from east to west, consistent with a 
known IR color gradient in the source. In NGC 3603, the H90a, He90a and the H113/3 
lines are detected from 13 individual sources. The Y + (He/H) ratios in the two sources 
range from 0.08±0.04 to 0.26±0.10. The H113/3/H90a ratio in NGC 3576 is close to 
the theoretical value, suggesting that local thermodynamic equilibrium (LTE) exists. 
This ratio is enhanced for most regions in NGC 3603; enhanced (3/a ratios in other 
sources have been attributed to high optical depth or stimulated emission. We compare 
the morphology and kinematics of the ionized gas at 3.4 cm with the distribution of 
stars, 10//m emission and H2O, OH, and CH3OH maser emission. These comparisons 
suggest that both NGC 3576 and NGC 3603 have undergone sequential star formation. 



Subject headings: HII regions: individual (NGC 3576, NGC 3603) — HII regions: 
kinematics and dynamics — ISM: abundances 



1. Introduction 



NGC 3576 (D = 3.0 ± 0.3 kpc) and NGC 3603 (D = 6.1 ± 0.6 kpc) are two of the highest 
luminosity optically visible HII regions in our Galaxy (Goss & Radhakrishnan, 1969). Indeed 
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NGC 3603 (with a bolometric luminosity of 10 7 L ) contains some 10 4 M Q of ionized gas, and 
aside from W49A, may be the most massive HII region in the Galaxy (Eisenhauer et al. 1998). 
Sources such as these are ideal laboratories to study the formation and evolution of massive stars 
and the initial mass function (IMF) of giant star forming regions. Similar regions that are visible 
with the Very Large Array (VLA) have been studied in detail. In particular, radio recombination 
line (RRL) studies of the Sgr B2 and W49A star forming regions with the VLA (Gaume et al. 
1995; De Pree et al. 1995: Mehringer et al. 1995) have revealed small-scale outflows, very broad 
recombination lines (AVfwhm ~80 km s _1 ), expanding shells of ionized gas, and extremely 
compact structures. NGC 3576 and NGC 3603 appear to be similar regions in terms of luminosity, 
total mass and gas kinematics. Existing optical studies of NGC 3603, for example, have revealed 
broad emission lines and a stellar wind bubble, indicative of the early stages of a star forming 
region coming into equilibrium with its environment (Clayton 1986). 

Though close in projection (~40'), the two regions are separated by about 30 km s _1 in radial 
velocity as indicated by HI absorption (Goss k, Radhakrishnan 1969) and hydrogen recombination 
line measurements (Wilson et al. 1970). Both regions are partially obscured by foreground dust, 
and studies of the kinematics of the associated ionized gas have been limited to optical wavelengths 
(Clayton 1986; Clayton 1990; Balick et al. 1980). Thus, existing velocity fields of the ionized 
gas are incomplete. The low resolution (4') 6 cm image of Goss &: Shaver (1970) shows that the 
radio emission associated with the two regions extends well beyond the detected Ha emission. 
In both sources, however, it appears that the large scale ionized structures may be driven by 
energetic processes occuring near their cores. Evidence from a variety of wavelengths (including 
the radio line and continuum data from this study) indicate that both regions may have undergone 
sequential star formation, in NGC 3576 from east to west, and in NGC 3603 from north to south. 

NGC 3603 and NGC 3576 have not been previously imaged with the Australia Telescope 
Compact Array (ATCA). Radio observations have been limited to low resolution line and 
continuum observations such as the ~l' observations of Retallack & Goss (1980) at 21 cm. 
Thus, our radio frequency observations {Ofwhm ~7") improve by factors of ~10 the best spatial 
resolution radio observations of the ionized gas associated with NGC 3576 and NGC 3603. 
For both NGC 3576 and NGC 3603, we present high resolution radio continuum images and 
radio recombination lines. New kinematic distances which are in reasonable agreement with the 
previously derived spectroscopic distances are derived from the observed line velocities and a 
revised Galactic radius of Rq=8.5 kpc. We examine the ionized gas morphology and kinematics in 
the two sources, and possible star formation histories of these two regions. In addition, we derive 
electron temperatures, H113/3/H90a ratios, helium abundances, and (where possible) carbon line 
strengths. 
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2. 



Observations and Data Reduction 



Radio continuum and recombination line observations were made for NGC 3576 and NGC 
3603 at the Australian Telescope Compact Array (ATCA) in the (750 m) B and (750 m) D 
configuration in two observing runs on the 13th and 19th of June in 1995 with a total of 24 hours 
of observing time. The observing frequency was 8.873 GHz, close to the the frequency of the 
H90a recombination line. The wide bandwidth (16 MHz) allowed the simultaneous observation of 
three additional recombination lines: He90a, C90a, and H113/3. Detailed observing parameters 
are presented in Table 1. The initial calibrations were made using the AIPS (Astronomical Image 
Processing System) software distributed by NRAO. Each data set was flagged and the flagged 
continuum data were then calibrated using the observed phase and flux density calibrators. These 
flags and calibrations were copied to the line data, where a bandpass calibration was applied and 
the continuum was subtracted from the u-v plane using the AIPS task UVLIN. Continuum images 
(with angular resolution of ~7") were made from the line free channels for each source using the 
AIPS task IMAGR. This spatial resolution corresponds to a linear size of approximately 0.09 pc 
for NGC 3576 (D= 3.0 kpc) and approximately 0.21 pc for the more distant NGC 3603 (D= 6.1 
kpc). The noise in the continuum images is dominated by the dynamic range limitations; thus the 
rms noise is significantly reduced in the continuum-subtracted line data. 

After the initial data reduction in AIPS, the line and continuum data were analysed using 
GIPSY (the Groningen Image Processing System). Integrated line profiles were made for each 
pixel above a continuum cutoff level for both NGC 3576 (5a = 90 mJy beam -1 ) and NGC 3603 
(5a = 55 mJy beam -1 ) using the GIPSY task PROFIL. Gaussian functions were then fitted to 
the spectral line data. The output from these fits are the line to continuum ratio (Tl/Tq), the 
integrated line strength, the central velocity (Vlsr) and the velocity full width at half maximum 
(AVfwhm) of the line emission in each region. The HII region NGC 3576 consists of a single 
extended source (containing several local bright spots) with an angular size of ~75" (~0.7 pc). In 
NGC 3603, thirteen individual sources were observed above a 5a level in the continuum image. 



Continuum images at 3.4 cm were made for both NGC 3576 and NGC 3603 using the AIPS 
task IMAGR. Figure 1 is the resulting 3.4 cm continuum image for NGC 3576. The crosses in 
the plot indicate the positions of the 10 /jm sources tabulated by Frogel & Persson (1974). Table 
2 lists the continuum parameters for NGC 3576. These parameters include peak flux density 
(Speak), integrated flux density (Stot), angular diameter in arcseconds (6), and linear radius in pc 
(r). Using the Parkes 64-m telescope, McGee k, Newton (1981) report a 14.7 GHz flux density of 
97.8 Jy in NGC 3576. We detect about 73% of the single dish flux density, or 71±7 Jy. From the 
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Results 



3.1. 



3.4 cm Continuum Emission 
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13 sources in NGC 3603, we detect a total flux density of 25 Jy, or 38% of the reported single dish 
flux density of 65.9 Jy. 

Figure 2 shows the 3.4 cm continuum image of NGC 3603. The crosses indicate the positions 
of the known 10 //m sources in NGC 3603 (Frogel et al. 1977). Each emission peak higher than 
the continuum 5a cutoff (55 mJy beam -1 ) has been designated as a separate region. The regions 
are labeled A to M, corresponding to increasing right ascension. The AIPS task JMFIT fitted 
two dimensional Gaussians to the individual continuum sources in NGC 3603. Table 2 presents 
continuum parameters for each region (A-M) in NGC 3603. These parameters include peak 
flux density (S pea k), integrated flux density (Stot), deconvolved geometrical mean diameter in 
arcseconds (6>), and the linear radius in pc (r). Table 3 presents derived physical parameters for the 
regions in NGC 3603. These calculated parameters include electron density (n e ), emission measure 
(EM), mass of ionized gas (Mhii), excitation parameter ([/), Lyman continuum photon flux 
(NLyc) and continuum optical depth (r c ). The continuum parameters were calculated assuming a 
uniform density, spherical, ionization bounded HII region, using the corrected formulae of Mezger 
& Henderson (1969). The calculated values are of course dependent on these assumptions, and in 
a complex region like NGC 3603, these assumptions may be overly simplistic. Figures 1 and 2 are 
plotted in B1950 coordinates. The last two figures in the paper (Figures 8 and 9) are plotted in 
J2000 coordinates to help in comparisons between the radio continuum emission and images at 
other frequencies, since images in the literature are shown in both epochs. 



3.2. Radio Recombination Lines 

3.2.1. NGC 3576 

A continuum-weighted, integrated line profile was generated for NGC 3576 using the GIPSY 
task PROFIL. The integrated profile was generated for all pixels above a continuum 5a cutoff 
level. In addition, integrated profiles were generated for each pixel above the continuum 5a level 
in a 7"x7" box centered on each of the positions of the Irs sources tabulated by Frogel & Persson 
(1974). The integrated source profile, and the four profiles line profiles were then fitted with 
Gaussian functions using the GIPSY task PROFIT. Figure 3a presents the integrated line profile 
and model fit for NGC 3576. Each of the four lines (H90a, He90a, C90a and H113/3) is clearly 
detected, and labeled in Figure 3a. The integrated spectra for the Irs sources (Irsl/2, Irs3, Irs4 
and Irs5; Frogel & Persson 1974) in NGC 3576 are shown in Figures 4a-d. The spectra in Fig. 3a 
and Figs. 4a-d show the data (points), Gaussian fits (solid line), and the residuals (dotted line) 
for each recombination line. The parameters of Gaussian fits to the lines (Vlsr, AVfwhm, and 
T;/T c ) are presented in Table 4. 

From the integrated emission from NGC 3576, we determine a helium abundance of 
Y + =0.09±0.01, and a H113/3/H90a ratio close to the theoretically predicted value of 0.276 
(/3/a=0.25±0.01). The electron temperature in NGC 3576 (T*) is calculated assuming local 
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thermodynamic equilibrium (LTE) from both the H90a and HI 13/3 recombination lines (using 
equation [22] from Roelfsema & Goss [1992]): 

T* = (6943z^ 1 ' 1 — - — — L_^0.87 K (!) 

with line width (AVfwhm) given in km s _1 and v in GHz. These LTE temperatures are given in 
Table 4. An electron temperature (T e ) may also be calculated with corrections made for non-LTE 
effects (using Eq. [23] from Roelfsema & Goss [1992]): 

T e = T:[b n (l-^)]°- 87 K, (2) 

where b n and [3 n values are from the tables of Salem & Brocklehurst (1979). 

For NGC 3576, a turbulent velocity (vturb) and an electron temperature (T e ) are calculated 
by solving two equations for line width (Eq. 3) using the hydrogen and helium (H90a and He90a) 
recombination line widths. The line width can be expressed as 

2 3AVl_3kT 1 
<v ^rb>~ 8ln2 M W 

where the vturb is the turbulent velocity in km s -1 , k is the Boltzmann constant, M is the mass 
of the atom, T e is the electron temperature, and AVjj is the line FWHM in km s _1 . Using this 
method for NGC 3576, we derive a turbulent velocity vt ur b=17±1.0 km s _1 and an electron 
temperature of (T e ) of 6100±200 K. 

In order to examine the spatial variation in the H90a line parameters, single Gaussian 
fits were made to each pixel above the continuum 5a level in NGC 3576. The results of these 
point-to-point fits to the H90a line are presented in Figures 5a and 5b. Figure 5a shows the 
spatial variation in the the central line velocity (Vlsr), with continuum contours overlaid on the 
greyscale representation of the ionized gas velocity. Figure 5b shows the spatial variation in the 
linewidth (AVfwhm) of the H90a line. For NGC 3576, a map of the LTE electron temperature 
was also generated, using Equation 1 (above) for each point above 10a in the continuum image. 
The resulting LTE temperature image is shown in Figure 5c, indicating a clear E-W temperature 
gradient. This temperature gradient has also been plotted using the AIPS task IRING, which 
averages emission in concentric semicircles. The semicircles (with a thickness of ~7"), are centered 
on the continuum peak of NGC 3576 and extend to the east. A plot of temperature as a function 
of distance from the peak (in arcsec) is presented as Figure 6. The apparent temperature gradient 
is discussed in §4.2.1. 



3.2.2. NGC 3603 

An integrated line profile was generated for the NGC 3603 region above a 5a cutoff in the 
continuum emission. Figure 3b presents the integrated line profile and model fit for NGC 3603. 
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Because of the large variation in line velocities across the source, the carbon line (apparent in 
some of the individual source profiles in NGC 3603) was fitted with the line width fixed at 10 
km s _1 . Gaussian fits were made to the integrated profile, and the parameters of these fits are 
given in Table 5. Figures 7a-7m present the radio recombination line profiles and Gaussian fits 
for each of the thirteen subregions (A-M) in NGC 3603. For all 13 subregions, the signal-to-noise 
ratio is sufficiently high to fit the H90a, He90a, and H113/3 lines. The C90a recombination line 
was successfully fit only in regions NGC 3603 B and H with a fixed carbon line width. Sources in 
which parameters were fixed and the fixed values are indicated parenthetically in Table 5. In a few 
cases, the fitted linewidth (AVfwhm) of the He90a line may be broad due to blending with the 
C90a line in the sources where a C90a line was not separately fitted. The line parameters from 
the Gaussian fits to the integrated emission from NGC 3603 and its 13 subregions are presented 
in Table 5. As with NGC 3576, the table also includes the LTE electron temperature (T*), the 
corrected electron temperature (T e ), the ionized helium to hydrogen ratio (Y + ), and the H113/3 
to H90a ratio (f3/a), and the turbulent velocity (vturb)- For NGC 3603, the non-LTE electron 
temperature was calculated using Equation 2 (above). The turbulent velocity was calculated from 
the width of the H90a line and the non-LTE electron temperature using Eq. 3 (above). These 
temperatures and turbulent velocities are given in Table 5. 

4. Discussion 

4.1. Distances to NGC 3576 and NGC 3603 

McGee & Newton (1981) detected recombination lines at 14.7 GHz from 25 high emission 
measure HII regions using the Parkes 64-m radio telescope. Both NGC 3576 (1109-610) and NGC 
3603 (1112-610) were included in their survey. McGee & Newton detected the H76a line from 
NGC 3576 and NGC 3603 at -24.1±0.1 km s" 1 and 7.7±0.1 km s" 1 respectively. In addition to 
H76a, McGee & Newton detected the He76a line toward both sources, and the carbon line toward 
NGC 3576. For NGC 3576, the detected line width in the H90a line (AV FWH m=28.7± 0.1) is in 
good agreement with the H76a line width published by McGee & Newton (AVpvF//M=29.1±0.1). 
The integrated H90a line width for NGC 3603 (AVfw.h\m=36.9±0.25) is in reasonable agreement 
with the value for the H76a line given by McGee & Newton (AVWhm=41.9±0.2). 

The distance to NGC 3576 has been discussed at length in two papers, Goss & Radhakrishnan 
(1969), and Persi et al (1994). Goss & Radhakrishnan found a kinematic distance of 3.6 kpc, just 
beyond the tangential point of the Galaxy, using the Schmidt galactic model and a value of Ro = 
10 kpc. This distance is derived from a line center velocity of -23 km s -1 , and circular motion 
within the galaxy. This kinematic distance agrees well with the findings from Caswell & Haynes 
(1987). Our new observations, using a revised Ro = 8.5 kpc, and a line velocity of -20.0 ± 0.1 
km s _1 suggest a distance of 3.0±0.3 kpc for NGC 3576. Because no ionizing cluster is optically 
visible, no reliable spectroscopic distance has been determined. Humphreys (1972), in a study of 
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the effect of streaming in the Carina arm, determined a distance of 3.2 kpc for NGC 3576, and 
detected both velocity variations (~4 km s _1 ) and a large spread in the individual velocities of 
early type stars in this direction. Persi et al. (1994) derive a distance based on the star HD 97499 
(position plotted in Figure 8). However, the high resolution radio continuum image of the source 
(Fig. 8) shows that there may not be a direct association between HD 97499 and NGC 3576. HD 
97499 (RA = 11 12 10.1, Dec = -61 18 45.1 J2000) is located -2' east of the peak in the 3.4 cm 
continuum emission. As a result, it is unlikely that HD 97499 contributes greatly to the ionization 
seen in NGC 3576. Finally, the stellar velocity of HD 97499 (-32 km s" 1 ± 11 km s' 1 ; Crampton, 
1972) does not agree well with the observed velocities in the ionized gas (-23 km s _1 [McGee &i 
Gardner, 1968] and -25 km s^ 1 [Caswell & Haynes, 1987]). Thus, distances associated with HD 
97499 should not simply be applied to NGC 3576. 

The determination of the distance to NGC 3603 is aided by the presence of an optically 
visible ionizing cluster. Distance estimates made from stellar observations can be compared to 
kinematic distances derived from the ionized gas. The earliest spectroscopic distance estimate was 
made by Sher (1965) at 3.5 kpc. Without a large number of faint stars visible, the spectroscopic 
distance estimate was admittedly crude. Shortly thereafter, Goss & Radhakrishnan (1969) found 
a distance of 8.4 kpc based on the observed line velocity of 10 km s _1 , using Ro = 10 kpc. Caswell 
& Haynes (1987) also found a kinematic distance using Ro = 10.5 kpc and V = +11 km s _1 of 
8.6 kpc. McGee Sz Newton (1982) report an H76a line velocity of 7.7±0.1 km s _1 . Our velocity 
for the integrated line in NGC 3603 is Vj,Sii=9.1±0.1, consistent with previously reported values. 
The H90a velocities found for the thirteen individual regions range between —2 km s _1 and 19 
km s _1 . Using the velocity from the integrated line (Vlsr=9-1 km s _1 ) and Ro = 8.5 kpc, we 
determine a kinematic distance of 6.1±0.6 kpc. 

UBV stellar photometry of NGC 3603 has resulted in distance estimates as low as 3.5 kpc 
(Sher, 1965) and as high as 8.1 kpc (Moffat, 1974). Melnick k Grosb^l (1982) derive a distance 
of 5.3 kpc. Moffat (1974), using the stars observed by Sher (1965) and a different reddening scale 
derive a distance of 8.1 ± 0.8 kpc. Finally, Melnick, Tapia & Terlevich (1989) use UB V photometry 
of the ionizing cluster of NGC 3603 to find a distance modulus of (m — M)o=14.3, which 
corresponds to a distance of ~7.2 kpc. Their value is in good agreement with both Moffat (1974), 
and the kinematic distances of Goss & Radhakrishnan (1969), and Caswell & Haynes (1987), 
and in reasonable agreement with the distance derived from our integrated radio recombination 
line velocity. The smaller distances inferred from the measurements of Sher (1965) and Melnick 
& Grosb(/>l (1982) are likely the result of errors in photometry due to nebular contamination 
(Melnick, Tapia & Terlevich 1989). 



4.2. NGC 3576 



NGC 3576 has been extensively observed in the radio continuum (Goss & Shaver 1970) and in 
radio recombination lines (McGee & Gardner 1968, McGee & Newton 1981, Wilson et al. 1970). 
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Recombination line studies of the region show the line velocities to be centered at approximately 
-23 km s _1 . The core (diameter ~ 1.5') of NGC 3576 is the location of five identified infrared 
(10/xm) sources (Frogel & Persson 1974), as well as CH 3 OH (Caswell et al. 1995) and H 2 
(Caswell et al. 1989) masers. The presence of bright 10/xm emission, water masers and compact 
thermal radio emission are typical indications of the early stages of star formation and dense 
circumstellar environments. 

Persi et al. (1994) have identified 135 individual stars with IR excess (to distinguish cluster 
stars from red field stars) in an extremely young galactic cluster associated with the HII region 
NGC 3576, using JHK images and photometry. Approximately 40 of these IR sources are 
members of the cluster. In their JHK images, Persi et al. (1994) detect a color gradient in the 
near-IR, indicating increasing extinction toward the western edge of NGC 3576. This gradient 
suggests that star formation may have progressed in NGC 3576 from east to west, with the 
youngest stars currently located at the western edge of the source. 

In addition, Persi et al. (1994) provide a summary of existing studies of the region, and an 
overlay of the 21 cm continuum (Retallack & Goss 1980), the 10 /xm continuum sources (Frogel 
& Persson 1974) and their K band mosaic of the region. The overlay (their Fig. 7) clearly shows 
that the radio continuum, near-IR continuum and 10/«n continuum are all strongly peaked toward 
the western edge of the source. Our 3.4 cm radio continuum image (Fig. 1) improves significantly 
on the resolution of Retallack k, Goss (1980), and shows that there is a coincidence between the 
location of the infrared sources (Irsl-Irs5) tabulated by Frogel & Persson (1974) and the extended 
ionized radio continuum emission. Figure 1 shows the 3.4 cm radio continuum with crosses 
and labels indicating the positions of the 10/mi peaks (Frogel et al. 1977). All five sources are 
coincident with the ionized gas emission and are near local peaks in the radio continuum emission. 

4-2.1. Characteristics of the Ionized Gas 

Our radio recombination line observations have allowed us to examine in detail the variation 
in electron temperature and helium abundance in NGC 3576. Our results in both cases seem to 
strengthen the case for sequential star formation in that source. The LTE electron temperatures 
show a clear gradient from east to west, with temperature increasing with distance (east) from 
the continuum peak. The gradient is apparent in Figures 5c and 6. Figure 5c shows the spatial 
variation in electron temperature across the source. Figure 6 is a plot of temperature as a function 
of distance from the 3.4 cm continuum peak. The LTE electron temperature increases with 
distance from 6000 K to a maximum of ~8,000 K, and then falls to lower levels in the more diffuse 
gas to the east. If the ionizing stars are those IR sources coincident with the radio continuum, 
this rising temperature signature is expected (Hjellming 1966). Lower energy photons on average 
would be absorbed first, with the higher energy photons penetrating deeper into the diffuse gas. 
Thus, the temperature gradient confirms a scenario in which the youngest, hottest stars in the HII 
region are located near the sharp western edge of the radio continuum. 
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Under LTE conditions, the theoretical value of the H113/3 to H90a hydrogen ratio is 0.276 
(Shaver & Wilson, 1979). Deviations from this theoretical ratio indicate local deviations from 
LTE, and the ratio can be used to correct the derived electron temperature. The HI 13/3 to H90a 
ratio in NGC 3576 varies from ~0.25±0.05 to ~0.35±0.05, with the lowest values associated 
with the bright peak in the radio continuum along the western edge of the HII region. Pressure 
broadening can lower the (3 /a ratio from the theoretical value. When the (5 /a ratio is lower than 
the theoretical value, electron temperatures derived from the (5 lines will of course be higher than 
those derived from a lines. In fact, for NGC 3576, the electron temperatures derived from the 
H113/3 line are ~11 % higher than those derived from the H90a recombination line. 

A small number of Galactic objects are known to have helium abundances in excess of 20%. 
Such elevated helium abundances have been reported in the W3 core (Adler, Wood & Goss 1996, 
Roelfsema, Goss & Mallik 1992), and are thought to be associated with the environments of 
evolved stellar objects that are rapidly losing mass. A small but significant positive gradient in 
helium abundance (Y + ) is detected from the western peak of the source (Y + = 0.08 ± 0.01) to 
the more diffuse ionized gas to the east (Y + = 0.14 ± 0.03). This range is consistent with the Y + 
value reported by McGee & Newton (1981) of N(He+)/N(H+)=0.091. The abundance gradient is 
consistent with the suggestion of Persi et al. (1994) who detect a color gradient in the near-IR 
indicating increasing extinction to the west of the source. Persi et al. (1994) propose that the 
youngest stars are located at the western edge of the ionized gas. The higher helium abundances 
found to the east of NGC 3576 may indicate the presence of a population of older stars and their 
associated mass loss. 

4.2.2. Gas Kinematics in NGC 3576 

Figures 5a and 5b show that both broad lines and velocity gradients are present in the ionized 
gas associated with NGC 3576. Double-peaked emission lines are located over a ~15" region at 
the northern edge of the source. In Fig. 5a, the double-peaked lines have been fitted with a single 
Gaussian with a AVfwhm = 50 km s _1 . Closer inspection (see Figure 1 inset) shows that locally 
the lines are possibly double-peaked. The double-peaked recombination lines in this region have a 
velocity separation of ~25 km s -1 , and may be due to either the expansion of an ionized shell of 
gas, or the outflow of ionized gas from a local source. The broad line region is located to the east 
of the Irsl/Irs2 peak (Frogel & Persson 1974). The region of the broadest radio recombination 
lines corresponds spatially to a distinct gap in the 10 fim emission between Irsl/Irs2 and Irs5. 

If the split lines are due to expansion, then the rate of expansion is ~13 km s -1 . Double-peaked 
lines with similar velocity separations have been seen in some of the compact HII regions in Sgr 
B2 North (De Pree et al. 1996). The spatial resolution of the current data is insufficient to 
determine whether an expanding shell of ionized gas is located at this position. Alternatively, the 
double-peaked lines may arise near the base of a large scale ionized outflow. 
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A large, approximately north-south velocity gradient of 30 km s _1 pc _1 is observed in the 
ionized gas, and is shown in Figure 5a. The velocity across the source varies from ~0 km s _1 
along the southern edge to ~-35 km s _1 in the north. The gradient is not likely due to rotation 
in the ionized gas, since these rotational velocities would imply an approximate enclosed mass of 
M=2xl0 5 Mq. The gradient is more likely due to an ionized outflow, and may in fact be the base 
of a much larger scale outflow. A wide field of view optical image of the region (Fig. 8; Persi et 
al. 1994) shows several large loops and filaments of ionized gas extending several minutes of arc 
to the north of NGC 3576. 

4-2.3. Sequential Star Formation in NGC 3576 

Figure 8 shows the positions of cluster members with infrared excess (crosses; Persi et al. 
1994) overlaid on the 3.4 cm radio continuum. Approximately 40 infrared excess sources have 
been identified. Eight of the stellar positions are coincident with the radio continuum emission, 
with the majority of the sources located to the north, east and south of the radio continuum. 
Only a single infrared excess source is located to the west of the radio continuum peak. This 
spatial distribution would be expected if star formation were proceeding from east to west into 
the molecular cloud. Water (H2O) and methanol (CH3OH) maser emission (plotted) are located 
coincident with the peak in the radio continuum emission. 

The location of molecular gas in NGC 3576 is also consistent with a triggered star formation 
scenario. The peak in the 1.0 mm dust continuum emission (Cheung et al. 1980) is located to the 
west of the bright radio continuum peak. The 1.0 mm emission traces the core of the molecular 
cloud located to the west of NGC 3576. Similar distributions of ionized and molecular gas are 
found near the Sgr B2 F HII regions (De Pree et al. 1996), G34.3 (Gaume et al. 1994), and others. 
In fact, the morphology of NGC 3576 is "cometary", and closely resembles a mirror image of the 
cometary source G34.3 (Gaume et al. 1994) in which a dense molecular cloud core is located at 
the eastern periphery of the ionized gas. 

Taken as a whole, the radio, infrared, optical and molecular imaging of NGC 3576 suggest 
that the oldest ionizing sources are located to the east of the region, and the youngest ionizing 
sources to the west. The distribution of stars, ionized gas, molecular gas, and water masers in 
NGC 3576 is reminiscent of simple scenarios of sequential star formation. In this picture, the star 
formation process slowly "eats away" at a dense cloud, so that the most evolved stars are located 
far from the cloud core. The ionized gas is found at the periphery of the cloud, and the youngest 
sources (protostars) are still embedded in the nascent molecular material, but reveal their presence 
with water maser emission. In many such regions in the Galaxy, the stellar population is not 
observable. The ability to observe stars, ionized gas and molecular gas in this distant source make 
it unique in the Galaxy. 
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4.3. NGC 3603 

NGC 3603 fills a gap in luminosity between a region like Orion, with ~1 O-type star, and 
R136 in 30 Doradus, with ~100 O-type stars (Eisenhauer et al. 1998). The ionizing core of NGC 
3603 (HD 97950) consists of approximately 20 O-type and Wolf-Rayet stars (Moffat, 1983). From 
photometric measurements, Melnick, Tapia <fe Terlevich (1989) have identified 181 sources (50 
with U, B, and V photometry) associated with the stellar cluster within NGC 3603. The authors 
found an average age of 2.5 ± 2 Myr and evidence that star formation propagated from north 
to south, ending at the bordering molecular cloud (Cheung 1980). Brandner et al. (1997) have 
identified a ring-shaped nebula and bipolar outflows from Sher 25, a cluster member of HD 97950. 
Sher 25 (a B1.5Iab supergiant) is located approximately ~20" north of the core of the cluster. 
Brandner et al. present evidence for association of this star with HD 97950 and compare the star 
to the progenitor of SN 1987A. 

The initial mass function (IMF) of HD 97950, the stellar cluster associated with NGC 3603, 
has been discussed in recent articles by Moffat et al. (1994), Hofmann et al. (1995), Zinnecker 
(1995) and Eisenhauer et al. (1998). Because of the proximity of NGC 3603, high resolution, 
deep studies allow detailed counts of stellar mass distribution to be made. Moffat et al. (1994) 
observed the central cluster with the HST/PC, studying its dense stellar core and the population 
of Wolf-Rayet stars. Moffat et al. (1994) conclude that in terms of density and stellar population, 
NGC 3603 is a 'Galactic clone' of R136 in 30 Doradus. Hofmann et al. derive an IMF power law 
of r =—1.4 ± 0.6 (£ ocM r ) for 30-60 M & which is similar to that of 30 Doradus. Using speckle 
masking techniques with the HST, Hofmann et al. (1995) found a similar slope for stellar masses 
of 15 to 50 M ofT =—1.6 ± 0.2. Both Eisenhauer et al. and Zinnecker use adaptive optics to 
extend the IMF to stellar masses less than 1 M . Some theoretical and observational studies have 
claimed that starburst systems are deficient in low mass stars, but Eisenhauer et al. (1998) find 
that the IMF follows a Salpeter power law with an index of < —0.73, and that there does not 
appear to be a low mass truncation to the IMF down to the observational limit of 0.2 Mq. From 
such properties of NGC 3603, it may be possible to extrapolate properties of larger extragalactic 
starbursts. 

Sequential star formation has also been suggested in NGC 3603, proceeding from north to 
south (Melnick, Tapia & Terlevich 1989). The presence of evolved objects at the north end of 
the source and embedded IR sources and a molecular cloud to the south qualitatively confirms 
a sequential star formation scenario (Hofmann et al. 1995). Investigations of NGC 3603 have 
proceeded on two fronts: radio and optical observations of the ionized gas, and optical and 
infrared observations of the exciting stars. While the stars and gas in the giant HII region are 
visually observable, these high resolution radio observations present the first unobscured view of 
the ionized gas on small scales. Melnick, Tapia & Terlevich (1989) have tabulated the positions of 
181 stars in NGC 3603 from their CCD UBV photometry. The positions of the 50 cluster member 
stars for which U, B, and V magnitudes are tabulated are shown in Figure 9. The figure shows 
the positions of these stars (crosses) plotted with the 3.4 cm radio continuum. Coordinates for 
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the stellar positions were derived from pixel coordinates using the known position of the B1.5Iab 
supergiant (Sher 25) as given by Brandner et al. (1997). 

Optical studies indicate both the presence of high velocity gas and of highly evolved stars. 
Clayton (1986) first investigated the velocity structure in the ionized gas located at the core of 
NGC 3603, and found both large scale motions (perhaps kinematic evidence of earlier supernovae) 
and two smaller shells, possibly wind-driven bubbles. The Ha and [Nil] line profiles were made 
along two lines across the core of NGC 3603. One possible explanation of the wind-driven bubbles 
was suggested by Moffat et al. (1994) who detected three Wolf-Rayet (WR) stars in their HST/PC 
observations of the cluster core. Hofmann et al. (1995) confirm the presence of three, perhaps 
four WR stars in the cluster core. Brandner et al. (1997) provide fascinating evidence that one 
of the sources believed to be an older member of the cluster (Sher 25, a B1.5Ia supergiant) is at 
the center of a clumpy ring of optical emission and a possible bipolar outflow. They compare the 
source to the progenitor of SN1987A. 

Caswell & Haynes (1987) detected strong hydrogen radio recombination line emission from 
NGC 3603. The source was also included in the recombination line survey of McGee & Newton 
(1981), and both hydrogen and helium (at 14.7 GHz) were detected. Our integrated H90a line 
width (AVpw r i/M=36.9±0.25) is in reasonable agreement with the value for the H76a line given 
by McGee & Newton (AVWhm=41.9±0.2). 

The present study of the ionized gas in NGC 3603 addresses a number of the outstanding 
issues in this source: 

1. Are there any regions of enhanced helium abundance, as might be expected in the presence 
of 3-4 WR stars? 

2. Is there any evidence in the ionized gas emission of the bipolar outflow associated with the 
B1.5 supergiant Sher 25 (Brandner et al. 1997)? 

3. How does the location of ionized gas compare to the known positions of infrared sources 
(Frogel et al. 1977), WR stars (Moffat et al. 1994 and references therein) and maser emission 
(Caswell et al. 1989, Caswell et al. 1998)? 

4-3.1. Helium Abundance and (3/a Ratios 

McGee & Newton (1981) derived a helium abundance in NGC 3603 of Y+ =0.069 (no 
uncertainty given). Our calculated helium abundance from the integrated line profile is 
significantly higher (Y + =0.13±0.01), but in their notes on individual sources, McGee & Newton 
(1981) indicate that their He76a line profile was rather noisy. 

Moffat et al. (1994) detect three prominent Wolf-Rayet stars in NGC 3603 in their narrowband 
Hell (4686 A) images. WR stars have been cited as the possible explanation of enhanced helium 
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abundances (Adler, Wood & Goss 1996). It is certainly reasonable to look for evidence of such 
enhanced abundances in the vicinity of known WR stars. In fact, the highest helium abundances 
in NGC 3603 are detected in source C, which has a value of Y + = 0.26 ± 0.10. Interpretation of 
this high value is complicated by the large uncertainty. The ionized gas in NGC 3603 is located on 
the periphery of the WR wind-blown bubble, so it is possible that ionized helium (Y + ) abundances 
are enhanced by the presence of an evolved star. 

The observed H113/3/H90a ratios for the 13 sources in NGC 3603 range from 0.29 ± 0.04 to 
0.53 ± 0.12, all above the theoretical value (/3/a=0.276). Most detections are within 2cr of this 
theoretical value, but a few sources (NGC 3603 E, H, I, and K) exceed the theoretical value by 
more than 3a. Anomalously high (3 /a ratios were found from H92a and H115/3 recombination 
line observations in two Galactic center regions (GO. 18-0.04; the Sickle and the Pistol; Lang, et 
al., 1997). Thum et al. (1995) propose a model that can explain enhanced (3 /a ratios as the result 
of high continuum optical depth. However, NGC 3603 A-M appear to have high (3/a ratios, and 
yet be low density (n e <10 4 cm~ 3 ) regions with r c <1. Cersosimo & Magnani (1990) suggest 
that under such conditions, high [3 /a ratios may result from stimulated emission caused by a 
background source. 



4-3.2. Ionized Gas Kinematics 

No ionized gas is detected above a 5a continuum level at the location of Sher 25. The optical 
study of Clayton (1986) indicates clearly that there are both large scale and small scale motions 
in the ionized gas. Plate 1 of Clayton shows an Ha image of the central 6' of NGC 3603, and the 
positions of the two slit positions through the core. The bright ionized gas observed in our high 
resolution image (see Fig. lb and Fig. 9) lies to the south of the cluster core. On scales probed 
by our high resolution observations, only one source (NGC 3603 F) has clearly double-peaked 
recombination line emission. Two Gaussians with velocities separated by ~29 km s~ 1 (Vlsr 
=27.5 km s" 1 and —1.1 km s _1 ) have been fitted to the H90a line in NGC 3603 F. The velocity 
separation is similar to that seen in other Galactic wind-driven bubbles and also in NGC 3576. 
In the radio continuum, source F is a relatively faint, marginally resolved source located near 
the western edge of the southern bright region. NGC 3603 F has a deconvolved diameter of 
&FWHM ~25", giving the source a linear diameter of 0.5 pc (D = 6.1 kpc). Thus, source F may be 
an expanding shell of ionized gas. A small optically observable wind-driven "stellar bubble" with 
r ~0.6 pc has been observed centered on the HD97950 star cluster in NGC 3603 (Balick, Boeshaar 
& Gull, 1980). 

In NGC 3603, the individual source velocities (as derived from the H90a line) range from 
— 3.7±0.3 km s _1 in source A to to 15.9±0.2 km s _1 in source I. The broadest lines are detected 
in NGC 3603 F (AVfwhm=52.7±2.2 km s _1 ). As indicated above, the broad lines in this source 
likely result from the double-peaked lines apparent in Fig. 7f. Broad lines can result from a variety 
of local conditions including pressure broadening or spatially unresolved velocity gradients. Given 
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the relatively high frequency of the observations and other evidence for velocity gradients in the 
source, it is likely that the broad lines in source M (located to the north of the cluster core) are 
related to gas kinematics. NGC 3603 M has tWft — 26 ± 10 km s _1 . Other localized regions in 
the source have very broad recombination line emission. Lines as broad as AV ~80 km s _1 are 
detected in the low surface brightness gas between sources NGC 3603 D and F. And in the bright 
core of NGC 3603 D (^0.07 pc in diameter) line widths approach AV ~75 km s _1 . In a sequential 
star formation scenario, sources D and F are located at the boundary of the molecular cloud 
and the stellar cluster, and would be in the early stages of their evolution, perhaps still coming 
into equilibrium with their environment or photoevaporating material from the molecular cloud 
interface. 

4-3.3. Evidence for sequential star formation in NGC 3603 

In Figure 2, the positions of ten 10/xm sources (Frogel et al. 1977) are plotted over the 3.4 cm 
continuum. Figure 9 shows the positions of cluster stars and known maser emission. Error bars 
for the absolute positions of the IR sources are ±10", and the size of the crosses in Fig. 2 indicates 
the absolute error in the positions of the infrared sources. In several cases (e.g. Irs 1, 2, 9, 14), the 
infrared sources are coincident with radio continuum sources. 

Perhaps the clearest evidence for sequential star formation in NGC 3603 is the relative 
location of the radio emission from the ionized gas, the stellar positions, and the position of 
embedded infrared sources and maser emission. The distribution is very much like that observed 
in NGC 3576. The most evolved stars are located to the north, centered in what appears at many 
frequencies to be the wind-blown bubble resulting from the WR stars located there. The molecular 
cloud core is located farthest to the south, and the ionized gas emission (apparent at optical and 
radio frequencies) is located at the periphery between the stars and the cold gas. 

5. Conclusions 

Both NGC 3576 and NGC 3603 provide rare opportunities to study massive star formation in 
a variety of spectral regions. Massive star forming regions like W49A and Sgr B2 are impressive at 
radio frequencies, but the stellar populations are inaccessible in the infrared. In these two sources, 
high resolution stellar photometry and radio continuum and recombination line observations 
provide a more complete picture of the physical conditions and complex gas kinematics. 

The following is a summary of our major conclusions: 

(1) In these high resolution radio continuum and recombination line observations, we have 
identified a number of previously unresolved radio continuum sources, particularly in NGC 3603. 

(2) The observed temperature gradient in NGC 3576 is in agreement with the known color 
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gradient in the source and supports a sequential star formation scenario that has been proposed 
for the region. 

(3) We have detected several compact sources in NGC 3603 with broad recombination lines. 
One source (NGC 3603 F) has double peaked profiles, with the peaks separated by ~30 km s _1 . 

(4) The brightest radio continuum emission in NGC 3603 is located to the south of the core 
of the stellar cluster HD 97950. There is no radio continuum or line emission apparent at the 
position of Sher 25. 

(5) Derived helium abundances and electron temperatures in NGC 3576 are in good agreement 
with values given in McGee & Newton (1981) while those in NGC 3603 are significantly higher. 

(6) Enhanced helium abundances ranging in value from Y + =0.16 to 0.26 are detected in 
NGC 3603 A-D and J. These regions may have had their helium abundance enriched by the 3 or 4 
known WR sources located at the cluster core. 

(7) In both sources, the relative positions of stars, ionized gas, molecular gas and masers is 
suggestive of a sequential star formation scenario. In NGC 3576 the star formation appears to 
have proceeded from east to west, in NGC 3603 from north to south. 
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Table 1: Observational Parameters of NGC 3576 and NGC 3603 



Parameter 



NGC 3576 



NGC 3603 



Date (configuration) 

Total observing time (hr) 

R.A. of field center (J2000) 

Decl. of field center (J2000) 

Major axis (") 

Minor axis (") 

Position angle (°) 

LSR central velocity (km s -1 ) 

Total bandwidth (MHz) 

Number of channels 

Channel Separation (kHz, km s -1 ). 
Spectral resolution (kHz, km s -1 )... 

Continuum noise (mJy beam -1 ) 

Line noise (mJy beam -1 ) 



13 June 1995 (750m D) 
19 June 1995 (750m B) 



13 June 1995 (750m D) 
19 June 1995 (750m B) 



12 


12 




n h 15 m 02 g. 


-61° 18' 31" 


-61° 15' 46" 


7.0" 


7.0" 


6.8" 


6.9" 


-9° 


-13° 


-23 


10 


16 


16 


256 


256 


62.5 (2.1) 


62.5 (2.1) 


75.6 (2.6) 


75.6 (2.6) 


18 


11 


2.1 


2.2 
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Table 2: Continuum Parameters for NGC 3576 and NGC 3603 



Source 




RA 






Dec 


Speak 


STot. 


9 


r 






(J2000) 




(J2000) 


(Jy beam -1 ) 


(Jy) 


(") 


(PC) 


NGC 3576 


U h 


ll m 51.32 s 


-61° 


18' 


43'.'5 


4.1±0.4 


71±7 


~90 


1.4 


NGC 3603 A 


ll h 


14 m 56.02 s 


-61° 


13' 


55'.'8 


0.10±0.01 


1.3±0.3 


24 


0.52 


B 




14 m 56.51 s 




13' 


17'.'0 


0.08±0.01 


2.5±0.4 


38 


0.80 


C 




15 m 00.90 s 




13' 


32'.' 7 


0.06±0.01 


0.85±0.3 


24 


0.52 


D 




15 m 01.99 s 




16' 


33'.' 


O.lliO.Ol 


1.2±0.3 


22 


0.46 


E 




15 m 02.62 s 




15' 


53'.' 8 


0.40±0.02 


1.9±0.2 


13 


0.28 


F 




15 m 06.13 s 




16' 


40'.' 5 


0.12±0.01 


1.7±0.3 


25 


0.52 


G 




15 m 08.76 s 




16' 


55'.' 7 


0.30±0.01 


4.4±0.3 


26 


0.54 


H 




15 m 09.45 s 




16' 


41'.' 4 


0.23±0.01 


4.6±0.3 


30 


0.64 


I 




15 m 10.36 s 




16' 


17'.' 4 


0.25±0.02 


1.3±0.2 


14 


0.30 


J 




15 m 14.24 s 




17' 


34'.' 


O.lliO.Ol 


2.2±0.3 


30 


0.65 


K 




15 m 18.75 s 




16' 


58'.' 1 


0.08±0.01 


1.2±0.3 


25 


0.53 


L 




15 m 24.20 s 




12' 


54'.' 


0.06±0.01 


0.94±0.3 


29 


0.61 


M 




15 m 31.59 s 




13' 


16'.' 7 


0.06±0.01 


1.1±0.3 


26 


0.55 
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Table 3: Derived Continuum Parameters for NGC 3603 



Source 


n e 

(10 3 cm" 3 ) 


EM 
(10 6 pc cm" 6 ) 


M H n 
(M ) 


U 

(pc cm -2 ) 


N LyC 
(10 48 s" 1 ) 


Tc 


A 


0.92 


0.9 


13 


50 


5.6 


0.005 


B 


0.67 


0.7 


35 


62 


11 


0.004 


C 


0.75 


0.6 


10 


44 


3.7 


0.005 


D 


1.0 


1.0 


10 


49 


5.1 


0.006 


E 


2.8 


4.6 


6 


57 


8.2 


0.027 


F 


1.0 


1.1 


15 


54 


7.1 


0.007 


G 


1.6 


2.8 


25 


76 


19 


0.016 


H 


1.3 


2.1 


32 


76 


20 


0.012 


I 


2.1 


2.7 


5 


50 


5.4 


0.016 


J 


0.87 


1.0 


23 


60 


9.5 


0.006 


K 


0.85 


0.8 


12 


48 


5.0 


0.005 


L 


0.66 


0.5 


15 


47 


4.8 


0.003 


M 


0.72 


0.6 


12 


45 


4.0 


0.003 
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Table 4: Recombination Line Parameters and Derived Quantities for or NGC 3576 





V 


AV 


Ti/T c 


T e * 




(km s" 1 ) 


(km s- 1 ) 




(K) 


H90a 


-20.0±0.1 


28.7±0.1 


0.105±0.001 


6300±700 


He90a 


-19.3±0.4 


24.8±0.8 


O.OlliO.001 




C90a 


-21.5±0.9 


8.75±1.8 


0.003±0.001 




H113/3 


-22.3±0.2 


33.6±0.4 


0.022±0.001 


7100±800 



H90a 

Irsl/Irs2 

Irs3 

Irs4 

Irs5 



-24.0±0.1 
-17.2±0.1 
-14.5±0.1 
-18.8±0.2 



29.1±0.1 
26.5±0.1 
32.0±0.3 
36.2±0.4 



O.lOOiO.001 
0.112±0.001 
0.081±0.001 
0.098±0.002 



6500±700 
6400±700 
7200±800 
5500±600 
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FIGURE CAPTIONS 

Fig. 1 The 3.4 cm continuum image of NGC 3576 is shown at a resolution of 7" (0.09 pc at 3.0 
kpc). The first positive and negative continuum contours are at 3a (54 mJy beam -1 ). Subsequent 
positive contours are 1.4, 2, 2.8, 4, 5.6, 8, 11.3, 16, 22.6, 32 and 45.3 times the 3a level. The 
peak continuum flux density is 4.1 Jy beam -1 The positions of the 10 //m sources (Frogel & 
Person 1974) are indicated with crosses. Errors in the absolute IR positions are ±4" (indicated by 
the cross size). The inset is an integrated profile made over a region of the source where broad 
recombination lines are detected. The best fit to the spectrum consists of two Gaussians separated 
by ~ 25 km s _1 . Epoch is B1950. 

Fig. 2 The 3.4 cm continuum image of NGC 3603 with the 13 named sources (A to M in order 
of increasing RA). The continuum data has a spatial resolution of 7" (0.21 pc at D = 6.1 ± 0.6 
kpc). The first contour is at 5a (55 mJy beam -1 ) with subsequent levels at 1.4, 2, 2.8, 4, 5.6, 8, 
11.3, 16, 22.6, 32 and 45.3 times the 5a level. The peak continuum flux density is 0.47 Jy beam -1 . 
The positions of the 10 /im sources (Frogel et al. 1977) are indicated with crosses. Errors in the 
absolute IR positions from this later study are ±10" (indicated by the cross size). Epoch is B1950. 

Fig. 3 (a) The continuum weighted, integrated line profile for NGC 3576. The recombination 
lines shown (left to right) are H113/3, C90a, He90a, and H90a. The solid line is the Gaussian fit, 
the crosses are the data points, and the dashed line is the residual. The spectral resolution is 2.6 
km s -1 . (b) The continuum weighted, integrated line profile for the bright southern region in NGC 
3603. The recombination lines shown (left to right) are H113/3, C90a, He90a, and H90a. The 
solid line is the Gaussian fit, the crosses are the data points, and the dashed line is the residual. 
The spectral resolution is 2.6 km s _1 . 

Fig. 4 The continuum weighted, integrated line profiles for the Irs sources in NGC 3576 (Frogel 
& Persson 1974), indicated in Figure 1. The solid line is the Gaussian fit, the crosses are the data 
points, and the dashed line is the residual. The spectral resolution is 2.6 km s _1 . Figures show (a) 
Irs 1/2, (b) Irs 3, (c) Irs 4, (d) Irs 5. 

Fig. 5 Line parameters of the H90a line in NGC 3576 overlaid on the 3.4 cm continuum contours. 
Fits were made to each pixel above the 5a level in the continuum, (a) The velocity (Vlsr) of 
the H90a line in NGC 3576 for each pixel is presented in grey-scale with the NGC 3576 3.4 cm 
continuum contours overlayed. The contour levels are as indicated in Figure 1. The grey scale 
velocity range is from -35 km s -1 to 10 km s _1 . (b) The full width at half maximum (AVfwhm) 
of the H90a line in NGC 3576 for each pixel is presented in false color. The line width range is 
from 20 km s -1 to 50 km s _1 . (c) The electron temperature calculated for each pixel in NGC 3576 
above 10cr in the continuum image. The temperature (grey scale) ranges from 5000 to 12000 K. 
The contours start at 4000 K and rise linearly to 12000 K in increments of 1000 K. 

Fig. 6 The integrated electron temperature in NGC 3576 for 18 concentric semicircles (width of 4 
pixels) centered on the peak in the 3.4 cm continuum image. Temperatures range from 6000 K to 
9000 K across the source. 
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Fig. 7 The integrated line profiles from the 13 individual HII regions in NGC 3603, as indicated 
in Figure 2. The recombination lines shown (left to right where detected) are HI 13/3, C90a, 
He90a, and H90<x The plot shows the Gaussian fit (solid line), data (crosses), and the residual 
(dashed line). The spectral resolution is 2.6 km s _1 . Figures 7a to 7m show the integrated profiles 
of sources NGC 3603 A-M respectively. 

Fig. 8 The 3.4 cm continuum contours for NGC 3576. The stellar positions of cluster member 
stars with infrared excess (Persi et al. 1994) are plotted as crosses. Also plotted are the positions 
of the star HD 97499 (diamond), and the H 2 (triangle; Caswell 1989) and CH 3 OH (box; Caswell 
et al. 1995) maser positions. Continuum contours are as indicated in Figure 1. Epoch is J2000. 

Fig. 9 The 3.4 cm continuum contours in the vicinity of the known stellar cluster HD 97950. 
The positions of ~50 cluster members with measured U, B, and V magnitudes (Melnick, Tapia & 
Terlevich 1989) are plotted as crosses. Coordinates of the 50 stars were derived from pixel offsets 
from the known position of Sher 25 (a B1.5Iab supergiant) as given by Brandner et al (1997). Also 
indicated is the position of HD 97950 (diamond), and H2O (triangles; Caswell 1989) and OH (box; 
Caswell 1998) maser positions. Continuum contours are as indicated in Figure 2. Epoch is J2000. 
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Table 5: Recombination Line Parameters and Derived Quantities for NGC 3603 



Source 




y 

(km s" 1 ) 


AV 

(km s- 1 ) 


Ti/T c 


T e * 
(K) 


(K) 


Y+ 


P/a 


v turh 


Total 


H90a 
He90o 
C90a 
H113/3 


9.14±0.1 
8.9±1.0 
6.2±1.9 
5.9±0.3 


36.9±0.25 
38.9±2.6 
10.0 (fixed) 
42.6±0.8 


0.075±0.001 
0.009±0.001 
0.002±0.001 
0.026±0.001 


6600±500 
4800±400 




0.13±0.01 


0.39±0.01 


24±2 


A 


H90a 

He90« 

H113/3 


8.08±0.2 
4.0±1.3 
2.5±0.8 


26.3±0.5 
25.0±3.1 
31.6±1.9 


0.107±0.003 
0.019±0.003 
0.033±0.003 


6300±1100 
4900±900 


6700±1200 
5200±1000 


0.17 ±0.03 


0.31±0.03 


14±3 


B 


H90a 
He90a 
C90a 
H113/3 


-3.65±0.3 
0.9±2.0 
-12.5±1.3 
-7.4±1.1 


28.3±0.7 
28.8±4.8 
10 (fixed) 
29.0±2.6 


0.103±0.004 
0.016±0.004 
0.015±0.003 
0.029±0.004 


6300±1600 
6000±1500 


6400±1600 
6100±1500 


0.16 ±0.04 


0.29±0.04 


17±4 


C 


H90a 

He90a 

H113/3 


3.8±0.7 
9.1±4.3 
-5.2±4.1 


33.2±1.7 
44.1±10 
49.1±10 


0.104±0.007 
0.020±0.006 
0.023±0.006 


4300±1700 
4400±1700 


4200±1700 
4300±1700 


0.26 ±0.10 


0.32±0.11 


22±9 


D 


H90a 

He90« 

H113/3 


-0.80±0.4 
1.1±2.2 
-1.7±1.0 


35.9±1.1 
32.8±5.6 
31.4±2.6 


0.078±0.003 
0.014±0.003 
0.029±0.003 


6400±2000 
5700±1800 


6800±2100 
6100±1900 


0.16 ±0.05 


0.33±0.05 


23±7 


E 


H90a 

Hc90a 

H113/3 


1.9±0.2 
-2.2±1.5 
-2.3±0.6 


41.9±0.5 
36.9 (fixed) 
46.5±1.4 


0.066±0.001 
0.008±0.001 
0.022±0.001 


6700±600 
5300±500 


8500±800 
6700±600 


0.11±0.01 


0.36±0.02 


27±3 


F 


H90a 

He90« 

H113/3 


11.5±0.8 
27.6±2.7 
5.9±2.4 


52.7±2.2 
22.7±6.9 
59.9±6.8 


0.056±0.003 
0.011±0.005 
0.020±0.003 


6500±3300 
4700±2400 


7000±3600 
5100±2600 


0.08±0.04 


0.40±0.08 


36±19 


G 


H90a 

He90o 

H113/3 


15.5±0.4 
20.0±1.6 
8.3±1.2 


25.6±1.0 
17.2±3.8 
31.2±3.1 


0.115±0.006 
0.024±0.001 
0.040±0.006 


6200±2200 
4400±1500 


7900±2800 
5600±1900 


0.14±0.05 


0.43±0.07 


13±5 


H 


H90a 
He90o 
C90a 
H113/3 


11.8±0.1 
12.4±0.9 
8.2±1.6 
8.5±0.4 


32.2±0.3 
29.3±2.2 
10.0 (fixed) 
37.9±0.9 


0.082±0.001 
O.OlliO.001 
0.003±0.001 
0.028±0.001 


6900±1200 
5100±900 


8300±1400 
6300±1100 


0.12±0.02 


0.40±0.02 


19±3 


I 


H90a 

He90o 

H113/3 


15.9±0.2 
7.9±3.1 
13.4±0.7 


34.8±0.6 
44.9±8.0 
39.6±1.8 


0.068±0.002 
0.006±0.001 
0.023±0.002 


7900±2200 
5900±1600 


9300±2600 
6800±1800 


0.11±0.03 


0.39±0.03 


21±6 


J 


H90a 

He90a 

H113/3 


7.3±0.3 
2.9±3.6 
5.3±1.2 


39.9±0.9 
59.0±9.7 
47.8±2.9 


0.085±0.003 
0.011±0.003 
0.031±0.003 


5700±1600 
3600±1000 


5800±1600 
3700±1000 


0.18±0.05 


0.43±0.05 


26±7 


K 


H90a 

He90o 

H113/3 


5.5±0.3 
7.9±1.9 
3.8±0.7 


40.8±0.7 
31.2±4.7 
43.8±1.9 


0.093±0.002 
0.012±0.003 
0.037±0.002 


5100±1000 
3600±700 


5000±1000 
3500±700 


0.10±0.02 


0.43±0.03 


28±6 


L 


H90a 

He90o 

H113/3 


6.4±0.6 
5.7±4.5 
2.7±1.3 


38.2±1.5 
35.0 (fixed) 
32.9±3.1 


0.142±0.008 
0.018±0.004 
0.063±0.008 


3800±1000 
2900±800 


3700±1000 
2800±800 


0.11±0.03 


0.37±0.07 


26±7 


M 


H90a 

He90o 

H113/3 


4.6±0.5 
8.5±6.3 
4.3±2.1 


37.3±2.1 
33.8 (fixed) 
44.9±5.4 


0.115±0.009 
0.015±0.005 
0.051±0.009 


3600±1300 
2600±1000 


3500±1300 
2500±1000 


0.11±0.04 


0.53±0.12 


26±10 




Velocity (km/s) 
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